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Abstract 

We perform a statistical analysis with the prospective results of future experiments 
on neutrino-less double beta decay, direct searches for neutrino mass (KATRIN) and 
cosmological observations. Realistic errors are used and the nuclear matrix element 
uncertainty for neutrino- less double beta decay is also taken into account. Three 
benchmark scenarios are introduced, corresponding to quasi-degenerate, inverse hi- 
erarchical neutrinos, and an intermediate case. We investigate to what extend these 
scenarios can be reconstructed. Furthermore, we check the compatibility of the sce- 
narios with the claimed evidence of neutrino-less double beta decay. 
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1 Introduction 



Neutrino mass and lepton mixing represent an unambiguous proof that the Standard Model 
(SM) of elementary particles is incomplete. Various experiments with solar pQ, atmospheric 
[2] and man-made [3l H] neutrino sources imply non-trivial lepton mixing angles, as well 
as non-zero and non-degenerate neutrino masses. Their values are extremely suppressed 
with respect to the masses of the other (electrically charged) fermions of the SM. The most 
prominent and often studied mechanism to explain the smallness of neutrino masses is the 
see-saw mechanism [5] . The neutrino mass scale is here inversely proportional to the scale 
of its origin. In addition, lepton number violation is predicted: neutrinos are Majorana 
particles. Searching for this property will be a crucial test of the see-saw mechanism, 
but also of other mechanisms leading to small Majorana neutrino masses. Possible phe- 
nomenological consequences of lepton number violation are the generation of the baryon 
asymmetry of the Universe [6] or, at low energies, neutrino-less double beta decay (Oz//3/3) 
[7j. This decay of certain nuclei, (A, Z) — > (A, Z + 2) + 2e~, which has not yet been ob- 
served, clearly violates lepton number by two units, and is intensively searched for [7] . We 
will assume here that light Majorana neutrinos are exchanged in the diagram responsible 
for Qvfift. In this case, the amplitude for this process is proportional to the coherent sum 

3 

m ee = X) rrii , (1) 
i=i 

where rrii are the individual neutrino masses and U is the leptonic mixing, or Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS), matrix. The absolute value of m ee is called the effective 
mass. The entries U e i can be written as U e \ = cos 8u cos 6*13, U e i = sin #12 cos 6*13 e ia and 
U e3 = sin 6*i3 e l/3 , where a and (3 are two currently unknown "Majorana phases" and 6* 12 ,i3 
are mixing angles. While 6 13 is constrained mainly by short-baseline reactor experiments, 
612 is probed by solar and long-baseline reactor neutrino experiments. Their current best- 
fit values as well as la and 3a ranges can be obtained from three-flavor fits, the result 
being [5] 

sin 2 6 12 = 0.32 (±0.02)1°;°* , sin 2 9 13 = 0+°- 019 > - 050 . (2) 

In what regards the neutrino masses, for a normal ordering one has m 3 > m,2 > Tn\ 
with m\ = m\ + AuIq and m\ = m\ + Am\. In case of an inverted ordering we have 
m,2 > rrii > m3 with m\ = m 3 + At?t,q + Am\ and m\ = m| + Am\. Here Am^ and 
Am\ are mass-squared differences with best-fit values and 3a ranges (7.9l£i)- 10- 5 eV 2 and 
(2.6i°)]6) ' 10~ 3 e V 2 , respectively [8]. Quasi-degenerate neutrino masses occur when m\ 2 3 ^> 
Am\,Am 2 & . If neutrinos are Majorana particles, all low energy neutrino phenomenology 
can be described by the neutrino mass matrix m v = U* m^ iag U* . It contains nine physical 
parameters. Seven out of the nine parameters of the neutrino mass matrix appear in 
\m ee \. Therefore, it contains a large amount of information, in particular if complementary 
measurements of some of the other parameters exist. We also note that all parameters 
of m v which do not influence neutrino oscillations show up in the effective mass. Those 
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are the the Majorana phases and, in particular, the individual neutrino masses (neutrino 
oscillations are only sensitive to mass-squared differences). For a review on the dependence 
of \m ee \ on the various neutrino parameters see refs. [TJ [9], [10] and references therein. In 
the present paper, in contrast to other works statistically analyzing future neutrino mass 
measurements including QisftP [HI EH H2J, [13J, [151 EE Ej, we focus on the neutrino mass 
scale, i.e. the value of the smallest neutrino mass. To this end we define three natural 
benchmark scenarios and investigate how future experiments may be able to constrain 
them. Our goal here is to combine as much mass-related information as possible. 



2 Observables related to neutrino mass 

Currently the strongest experimental limits^ on the half-life of neutrino-less double beta 
decay are (all at 90 % C.L.) 1.9 • 10 25 y for 76 Ge [19] (see also [20]), T 1/2 > 3.0 • 10 24 y for 
130 Te [21], T 1/2 > 5.8 • 10 23 y for 100 Mo and T 1/2 > 2.1 • 10 23 y for 82 Se [22]. The existing 
limits on Ti/ 2 will be improved considerably (by two orders of magnitude or more) in the 
near future by various experiments [7] . The uncertainty in nuclear matrix element (NME) 
calculations is a serious problem to translate these bounds into upper limits on the effective 
mass [23l \Y7\ . We will take into account in particular this uncertainty in our analysis. 
Depending on the nuclei and NME, the current limit on the effective mass as extracted 
from the half-lifes given above lies between several tenths of and a few eV. This has to 
be compared with the predictions which can be made for the effective mass. Inserting the 
known ranges of the oscillation parameters, and varying the unknown parameters within 
their allowed ranges, one can generate plots as the ones in fig. [H 

They display (for U e 3 = 0) the effective mass as a function of the smallest neutrino 
mass, the sum of neutrino masses 

3 

s = H m i ( 3 ) 
i=i 

and the kinematic neutrino mass 



El^l 2 ^ 2 - (4) 

The latter two quantities can be measured through cosmological observations [21] and 
experiments like KATRIN |25J, respectively. The latter experiment has a 5<r discovery 
potential of 0.35 eV for nip, and a null result will lead to a 90 % C.L. limit of 0.2 or 0.17 
eV |26j . In the sensitivity range of KATRIN, the relation "imp = £ holds to a very good 
precision. Cosmology is expected to probe values of £ down to the 0.1 eV range [21] (to be 
specific, we take a value of 0.15 eV in fig. [[]). To achieve such impressive results, one takes 

1 We note that there is a claimed positive signal for 0v(3f3 from ref. [IS]. We will turn to this issue later 

on. 
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Figure 1: The effective mass as a function of the smallest neutrino mass, the sum of neutrino 
masses £ and the kinematic neutrino mass mp. The value = and the current 3a 
ranges of the other oscillation parameters have been used. 



advantage of future observations of weak gravitational lensing of galaxies, and the cosmic 
microwave background or detailed analyses of the 21 cm hydrogen emission lines at high 
redshift. It is fair to say that a conservative limit on E is 1 eV. This value corresponds 
roughly to the bound obtained from WMAP 5- year data alone [28]. Recall that neutrino 
mass bounds from cosmology depend strongly on the data sets, the priors and the model, 
i.e., adding parameters which are degenerate with neutrino masses will relax the bounds, 
see, e.g., [29] . Finally, current limits for trip are 2.3 eV [27]. 

The blue and yellow bands in fig. [I] correspond to the normal and inverted mass or- 
dering of the neutrinos, respectively. The darker areas in the blue and yellow bands are 
obtained when the oscillation parameters are fixed to their best-fit values and only the Ma- 
jorana phases are varied. The lighter areas correspond to the 3a ranges of the oscillation 
parameters. Note that this broadening is very weak for the maximum value of \m ee \ in 
the case of inverted mass ordering and for quasi-degenerate neutrinos. This is because the 
upper limits on \m ee \ are roughly \J Am\ and m^, respectively, and varying the oscillation 
parameters has very little impact. In the first plot of fig. [H, we have indicated three special 
values of |m ee | which correspond to the goals of the three phases of the GERDA experiment 
(where a certain NME has been assumed, see [30] for details). 
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3 Statistical analysis 



Now we will perform a statistical analysis to investigate how well it will be possible to re- 
construct different realistic physical scenarios with upcoming neutrino mass experiments. 
Note that, since we want to investigate realistic situations, we concentrate only on cases 
that can be probed in the near future. For definiteness, we consider the inverted mass 
ordering and three different scenarios called QV (quasi-degenerate), XJ\fT (intermediate) 
and X7i (inverted hierarchy) that are defined by different values of the smallest neutrino 
mass 7713. Note that the QV scenario would, to very large extent, also apply to a nor- 
mal mass ordering. The hypothetical "true values" for the different observables in these 
scenarios are: 



Scenario 


m 3 [eV] 


\m ee \ 


[eV] 


[eV] 


£ [eV] 


QV 


0.3 


0.11 - 


0.30 


0.30 


0.91 


ZA/T 


O.f 


0.04- 


O.ff 


(o.n) 


0.32 


in 


0.003 


0.02- 


0.05 


(0.05) 


(O.fO) 



We have used here the best-fit values for the oscillation parameters. The range for \m ee \ 
originates from the variation of the Majorana phases a and f3. Note that the KATR1N 
experiment will only be able to measure nip in the case of the QV scenario, while for I NT 
and X7i it will only provide an upper limit. The same is true for the measurement of E in 
the XTi scenario. These cases are indicated in the table by writing the respective values in 
brackets. 

Let us now give a summary of the different experimental errors and theoretical uncer- 
tainties. Regarding the error on the effective mass in 0v(3/3, we have to distinguish between 
experimental and "theoretical uncertainties", where the latter result from the NME uncer- 
tainty. The experimental error can be included by noting that the decay width depends 
quadratically on the effective mass. Thus, 

^(|m«U) = ^^^*!), (5) 

" J- obs 

where |m ee | ex p is the measured value of the effective neutrino mass and 0"(r oD s) is the 
experimental error on the measured decay width r obs for neutrino-less double beta decay. 
For definiteness, we choose the ratio of the latter two as 

a(r ° bs) -23.3%, (6) 



r, 



obs 



which is the value obtainable in the GERDA experiment [30]. We combine, similarly to 
the procedure developed in ref. the experimental error with the theoretical NME error 
via 

a(\m ee \) = (! + £) (jm ee | + o-(\m ee \ exp )J - \m ee \ , (7) 
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where £ > parameterizes the NME uncertainty and c"(l m ee| exp ) is given in eq. (151) . Fol- 
lowing ref. |14j . we define a covariance matrix 



c r 2/ \ , ^6 2 

Sot = <U ex (a) + 2^ ~ ~ a. 



dxj dxj 



where T% = \m ee \, T2 = X and T3 = m|. Furthermore, <r 2 (a) is the error on T a , and a, 6 
label the entries in the covariance matrix. The Xi are the oscillation parameters that enter 
\m ee \ (and mp, though in the observable range of mp they have basically no influence). The 
errors on the T a are given by eq. ((7|) as well as by a(m 2 ^) = 0.025 eV 2 [25| [26] and 
ct(£) = 0.05 eV [23]. 

Defining t> a = T a — (T a ) exp , where (T a ) ex p denotes the experimental value of T a , our 
X 2 -function to be minimized is 

X 2 = v T S~ l v. (9) 

All oscillation parameters are set to their current best-fit values and their (symmetrized) 
standard deviations are determined from their la-ranges, which is a good approximation 
for future 3<r-ranges. Anyway, the impact of different numerical values here would not lead 
to qualitatively different results. We first minimize the \ 2 from eq. (Q with respect to the 
Majorana phases a and (3. The resulting function is x 2 es = min ajj a % 2 . We then continue by 
plotting the resulting la, 2a and 3a ranges for the smallest neutrino mass 771,3 determined 
by setting A% 2 = x 2 es — x 2 esmin equal to 1, 4 and 9. This corresponds to a x 2 -f unc ti° n 
with one free parameter (namely m 3 ). |m ee | oxp is the assumed measured value of \m ee \, on 
which the reconstructed range of m 3 depends. The minimum in the \m ee \ exp -m 3 plane is 
determined such that A% 2 is zero in the true region of the corresponding scenario (e.g., 

OP). 

The results of our analysis are shown as the solid lines in the left column of figs. El 
[3] and HI In all cases, we have calculated the result for a consistent measurement (i.e., 
mp and S are measured at their true values in the corresponding scenarios). The NME 
uncertainties we have chosen are Q = (no uncertainty), 0.25 and 0.5. We have checked 
that values of ( > 0.5 will lead to results not too much different from the ones for £ = 0.5. 
The value £ = 0.25 is a quite typical one, cf. refs. [231 HZ]. This uncertainty arises from 
the highly non-trivial calculations of the nuclear part of the neutrino-less double beta 
decayprocess. Different methods, and even different Ansatze within the same framework, 
differ in their result, and their spread is commonly taken into account as "theoretical 
uncertainty". Glancing at Fig. 5 in ref. [31], where the results of different methods of the 
NME calculation are compared for different nuclei including 76 Ge, one can indeed see that 
the spread of the respective values around their mean value is about 0.2. We conclude that 
the values we use are realistic and typical. 

The true values of \m ee \ and 777,3 are marked by the vertical black lines. The plots 
illustrate how well we can reconstruct the different scenarios for the various values of NME 
uncertainty. Having a look at fig. [21 we see that the QT> scenario can be reconstructed 
quite well, which is not surprising since in that case the KATRIN experiment as well as the 
cosmological measurement will provide a non-trivial signal. E.g., for |7n ee | exp = 0.20 eV, 
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im [cV] m [oV] m 3 [eV] 



Figure 2: la, 2a and 3a regions in the m^-\m ee \ plane for the QV scenario. The 
left column shows the correct (solid line) as well as two possible incorrect cosmological 
measurements (dashed lines). The less desirable case, namely only taking into account a 
KATRIN measurement, is shown in the plots on the right. The area denoted HDM is the 
range of \m ee \ from the claim of part of the Heidelberg- Moscow collaboration. 



the 1, 2 and 3a ranges for m 3 are 0.28 - 0.32 eV, 0.27 - 0.33 eV and 0.25 - 0.35 eV, while 
the true value is 0.30 eV. Therefore, the reconstruction is quite accurate. This remains 
true also if the uncertainty in NME is non-zero because the plots are still narrow around 
the true value of m 3 (the numerical values suffer nearly no change) even though, with a 
larger NME uncertainty, also higher values of |m ee | exp are plausible. This is true for all 
three scenarios under consideration. 

Similar statements hold for the XMT scenario shown in fig. [31 even though mp cannot be 
measured now. However, because there will still be a measurement of E, we have sufficient 
information on the neutrino mass. In case the central measured value is \m PP \ = 0.08 eV 

i cc i exp 

and C = the ranges are 0.08 — 0.12 eV at la and 0.05 — 0.15 eV at 3a. In case of ( = 0.5 
we find 0.08 — 0.12 eV at la and 0.04 — 0.15 eV at 3a. The mass scale has now a 3a 
uncertainty of 50 %, to be compared with roughly 15 % in the QT> scenario. 

For 27i, in turn, there is no measurement that gives information m^. Hence, it is only 
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Figure 3: Same as fig. [2] for the 1MT scenario. 



possible to give an upper limit on the smallest neutrino mass, as illustrated by the long 
horizontal band in the left column of fig. HI Note that this band corresponds to the yellow 
band marking the inverted mass ordering in the upper plot of fig. [U This upper limit is 
almost trivial, i.e., it corresponds to the neutrino mass limit obtainable from Ovflft alone. 
To give some numerical values, for |<Ti e e| exp = 0.04 eV one would have the 1 (3)cr ranges 
m 3 < 0.03 (0.07) eV for ( = and for £ = 0.5. Due to the bound on E, there is very little 
dependence on £. 

Up to now, the discussion has focused on the case in which all measurements are 
compatible. As an example for inconsistency we discuss here a possible clash between 
results from KATRIN and from cosmology. To this end we leave (mg) exp equal to the true 
value of the corresponding scenario (new physics is not expected to influence rrip [32J) and 
take values of E exp which are smaller or larger than the true value. There are many scenarios 
or models in the literature which can lead to wrong values of E, see, e.g., refs. [33J. The 
result is shown by the areas within the dashed lines in the left columns of figs. [20 Having a 
look at QD first, we realize immediately that the physical range is reconstructed incorrectly. 
Hence, if there are systematic errors in the cosmo logical measurement, or unknown features 
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Figure 4: Same as fig. [2] for the X?i scenario. 



in cosmology which we are not aware of, a wrong neutrino mass is reconstructed. In the 
QT> case there is still information from KATRIN, which leads to a reconstructed neutrino 
mass at most one order away from the true value, even if the wrong E is taken into account. 
For the TNT scenario, however, there is no information from KATRIN. Consequently, it 
might be that a wrong upper limit on m 3 is concluded, as illustrated by the long band for 
^exp — 0.05 eV in the upper left plot of fig. [3j This is an example wherein one could draw 
a wrong conclusion by taking the cosmological measurement at face value. As expected, 
even worse cases may exist for the T7i scenario. E.g., in the upper left plot of fig. H]one 
would, for E exp = 0.3 eV, reconstruct a smallest neutrino mass of roughly 0.1 eV, to be 
compared with the true value of rris = 0.003 eV. For the TH, scenario, one might not even 
realize that there is an inconsistency, since in that case, the KATRIN experiment can only 
provide an upper limit which is too far away from the true value of m^. 

One possible cross-check (or the possible consequence if one indeed finds that the results 
from KATRIN and from cosmology do not fit together) would be to dismiss the cosmological 
data altogether. We have also analyzed this case. Here, S a b from eq. (jSJ) as well as v a 
would change from 3-dimensional to 2-dimensional objects while the rest of the procedure 
remains the same. The results for this analysis are plotted in the right columns of figs. [2]|H 
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again for different values of the NME uncertainty. For QV, the most optimal scenario, 
neglecting cosmology, would simply increase the errors in the determination of m 3 : e.g., 
for \m ee \ exp = 0.20 eV and C = the ranges are 0.26 - 0.34 eV at la and 0.16 - 0.41 eV 
at 3a, while for ( = 0.5 we find 0.26 - 0.34 eV at la and 0.13 - 0.41 eV at 3a. The NME 
uncertainty has now a slightly bigger impact, and the error on 7773 increases by a factor 
of three, since now it is about 50 % while it was roughly 15 % when E has been included 
in the analysis. For the TNT scenario, however, there is a major difference to the former 
case: since now there is no other measurement besides |m ee | exp providing information on 
7?i,3, we can only derive an upper limit instead of determining a certain range for 7773. This is 
indicated by the band in the upper right plot of fig. [3j Finally, for Z7i, the limit on m 3 gets 
only slightly worse compared to the case of a E, which is too small to be measured. In this 
case there would not even be a real drawback in taking into account the KATRIN result 
only. It remains to be said that in all cases a higher uncertainty for the NME does not 
significantly modify the conclusions in what concerns the value of 7773. Finally, it is worth 
mentioning that if in QT> scenarios the error on E is decreased (increased), the obtained 
error on the neutrino mass is decreased (increased) by approximately the same factor. 

With our analysis we can also compare the compatibility of our three benchmark sce- 
narios with the range for m 3 of 0.15 — 0.46 eV, calculated as the (global fit) 2a range in 
ref. [16] from the claim in ref. [18]. We give the implied range for m 3 as the gray band in 
figs. El [3] and HI We see that scenario QT> is consistent with the claim, even for a measure- 
ment of E = 0.6 eV, to be compared with the true value E = 0.9 eV. The TNT scenario 
(T7i scenario) is barely (very) incompatible for measured "true" values, but a too high 
value of E exp can lead again to compatibility. We see that testing the claim and comparing 
it with cosmology is a non-trivial task (see also [T7]). 

4 Conclusions 

In this work we have investigated possible constraints on the neutrino mass in future 
experiments. We assumed realistic errors on the observables, in particular for neutrino- 
less double beta decay. Then, we have checked how certain realistic benchmark scenarios, 
which correspond to different regimes for the smallest neutrino mass, can be reconstructed 
from future measurements. Furthermore, we have pointed out how wrong conclusions could 
be drawn from inconsistent results, i.e., if cosmology provides a wrong value for the sum 
of neutrino masses. In case of consistent measurements we may summarize as follows: 
typical 3a errors for quasi-degenerate neutrino masses range from roughly 15 % (including 
E) to 50 % (excluding E), where NME uncertainties play a larger role in the latter case. 
Intermediate scale masses can also be determined with 50 % uncertainty. In case of an 
inverted hierarchy, the effective mass is constant for a large range of the smallest mass, 
which allows only to derive upper limits on it. 
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